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Abstract 

At high temperatures the Aq component of the Yang-Mills field plays the role of the 
Higgs field, and the 1-loop potential V(Aq) plays the role of the Higgs potential. We 
find a new stable vortex solution of the Abrikosov-Nielsen-Olesen type, and discuss 
its properties and possible implications. 



Recently there has been renewed interest in the quantized Z(N C ) vortices as possible 
candidates for the confinement mechanism |IL EL H|. If vortices are physical objects (and 
not lattice or gauge artifacts) playing a role in the dynamics of the vacuum fluctuations 
they have to be found as stable solutions of the effective action obtained from integrating 
out high frequencies of the field. In the zero-temperature case the one-loop effective action 
for vortices has been introduced in our previous work [|J. The zero-derivative term of the 
effective action, i.e. the 'potential energy' of the vortex has been found in that paper: it 
does not lead to any stable solution. This is probably not surprising since all derivatives 
of the effective action should be summed up before one gets to a definite conclusion on the 
existence of vortex solutions. 

In this letter we address a related but simpler question on whether there are stable 
vortex solutions in the case of nonzero temperatures. Moreover, we restrict ourselves to 



a simple and more 'pure' case of 2 + 1 dimensions. If the temperature T is high enough 
vortices can be viewed as cylinders pointing in the 'time' direction, with a nontrivial profile 
in the transverse 'spatial' plane. We shall show that the Abrikosov-Nielsen-Olesen-type 
vortex solution indeed exists in this case, and we shall present its profile and energy per unit 
length, as functions of temperature. 

The 2 + 1-dimensional YM theory has a coupling constant g\ of the dimension of mass; 
high temperatures means T 3> g\ . One can always choose a gauge with the 'time' component 
of the YM field A^(x) independent of 'time'. In this gauge, one can integrate out the high- 
momenta components of the fields to obtain the low-momenta effective action for time- 
independent fields A*(x), (jl = 0, 1,2. The separation scale between high and low momenta 
is given, naturally, by the temperature: high momenta means p > 2nT. The effective action 
is obtained by integrating over nonzero Matsubara frequencies, u n = 2irnT, n > 0, and 
over zero (n = 0) Matsubara frequency but large spatial momenta. It can be expanded 
in powers of the spatial derivatives of the A a ^(x) fields, divided by appropriate powers of 
T. The zero-derivative term is the 'potential energy' V(Aq). Since the spatial size of the 
vortex solution is expected to be larger than 1/T one can neglect all terms of the derivative 
expansion except the first, zero-derivative term V(Aq). The latter can be written as || Q 



V(A ) = — sin 2 ixv [°°dp - P 2 °° sh P = rT 3 u 2 ( In 1 + 2.67575") + 0(z/ 4 ) 
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It is a periodic function of the dimensionless variable v with unit period, depicted in Fig.l. 
Notice that the potential energy (0) is nonanalytic at the minima v = integer. The fact 
that the second derivative of V{Aq) has a logarithmic singularity at v = integer is related to 
the infrared divergency of the Debye mass in 2 + 1 dimensions || [7j, and will have important 
consequences for the vortex solution. 

Adding V(A ) to the classical YM action, F" /%!, we get the energy functional in 
the spatial transverse plane: 
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where i,j = 1,2 denote spatial components, Di is the covariant derivative, Df 1 = di5 ab + 

Eq. (fj) presents a 2-dimensional YM + Higgs system where Aq plays the role of the 
adjoint Higgs field, with a peculiar form of the potential energy (|1]). 

We choose the following vortex-type Ansatz for the fields, restricting ourselves to the 
SU (2) colour group for simplicity: 
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Potential V(A_0) in 2+1 dimensions 
at nonzero temperatures 
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Figure 1: 2 + 1- dimensional periodic potential V(A ), eq. 

where p(p) and u(p) are trial profile functions of the distance p from the vortex center. 
The Ansatz (|3]) corresponds to the radial component A a p = and the azimuthal component 
A%{p) = 5 a3 p(p)/p. The closed Wilson loop in the representation labelled by spin J, circling 
around the center of the vortex in the transverse plane at distance p, is 
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where p(p) is the magnetic field flux along the vortex inside a tube of radius p. If p(p) — > 1 
at large p the Wilson loop Wj(p) — > (— 1) 2J . In particular, in the fundamental representation 
one has W1/2 — > — 1. This is the definition of the quantized Z(2) vortex. Our solution will 
be precisely of this type. 

Let us introduce the dimensionless radius 



x = P9s 
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and the dimensionless parameter 
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At high temperatures a < 1. In terms of these quantities the action functional for the 
vortex of length 1/T becomes 
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where we have introduced the dimensionless 'Higgs potential' 
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The Euler-Lagrange equations of motion for the profile functions p, z/(p) are: 
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We look for the solution with the following boundary conditions: 

fi(0) = i/(0) = 0, /i(oo) = i/(oo) = 1, (13) 

corresponding to the quantized Z(2) vortex, with the 'Higgs field' A going from the trivial 
minimum at p = to a non-trivial one at p — ► oo. 

The behaviour of the profile functions near the origin and at infinity can be found ana- 
lytically. At small x we get: 

p(x) = CiX 2 + c 2 x A + ... , 

v{x) = d\x + d 2 x z Inx + d$x 3 + ... . (14) 



The coefficients c±, d\ are arbitrary but the higher coefficients are determined from eqs. ( ID , 
0): 

c 2 = L ; d 2 = - ( ^ 1 , d 3 = —[-2c 1 + an(3 + 2h-4\nd 1 )], h= 1.67575... (15) 
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At large values of x the analytical solution is: 



p(x) ~ 1 — eiv^rexp ( — 7= ] , (16) 



a 



v(x) ~ 1 — fi exp (— Ttax 2 ^ , (17) 
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Profiles of Z(2) vortices 
in high-temperature 2+1 dimensions 
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Figure 2: Profile functions of the vortex for two values of temperature, corresponding to 
a =0.15 and 0.5. 



where the constants ei, f\ are not determined by the equations. Notice that the 'Higgs field' 
v(x) approaches its asymptotic value at infinity not as an exponent but as a gaussian. This 
is a consequence of the logarithmic divergence of the Debye mass in 2 + 1 dimensions. Being 
rewritten in original notations eq. (O) reads 
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We notice that the characteristic scale for the variation of A (p), namely l/g^y/T, is, at high 
temperatures, much larger than 1/T. It justifies neglecting derivative terms of Aq in the 
effective action and leaving only the zero-derivative term V(A Q ), as in eq. (0). 
The coefficients Ci,d 1 ,e±, fi are determined by solving eqs. ( fLOl , 



TT|) starting from the 



expansion ([14]) towards larger values of x and starting from the asymptotic form (|i6|J17l) 
towards smaller values of x, and matching the functions and their derivatives at some inter- 
mediate point x ~ 1; this is done numerically. The resulting profile functions p(x) and v(x) 
are shown in Fig. 2, for two values of a = 0.15 and 0.5. 

Integrating the profile functions over the whole range of x we get the action of the vortex 
SVortex = E±/T plotted in Fig.3, as function of temperature T. At T > gl the action is 
large, so that the vortices are exponentially suppressed. This is a theoretically 'clean' case: 
as explained above, the use of the energy functional (Q) for finding the vortex solution is 
justified. The vortices are in fact short and broad cylinders oriented in the 'time' direction. 
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Vortex action 
as function of temperature 
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Figure 3: The vortex action as function of dimensionless temperature T j g\. 



Indeed, the length of the cylinders is 1/T while the characteristic radius where the 'Higgs 
field' A reaches its asymptotic value 2nT is of the order of l/g$y/T ^ 1/T, see eq. (I7p. It 
should be mentioned, however, that the magnetic flux p(p) reaches its asymptotic value of 
1 at a smaller radius p ~ 1/T, see eq. (|16|). 

To build the 2 + 1-dimensional vacuum at high temperatures out of the vortices, one 
needs first of all to gauge rotate the 'Higgs field' Aq @) to, say, the third direction in colour 
space: one cannot add up vortices with different colour orientations of Aq at infinity. This 
is achieved with a help of a ^-dependent gauge transformation 



a; = rfWAjJM + is^-u^^u^), u{4>) = ( ^ X-^ 2 ) • (19) 

Under this gauge transformation the Aq component of the YM field becomes proportional to 
r 3 , A'q = 5 a3 2nT u(p). The gauge transformation (|19D is, however, discontinuous at = 2n, 
therefore the azimuthal component A', will now have a 'Dirac-surface' singularity at = 0: 

A'; = 5 a2 -[l-p(p)-27r5(<t ) )}. (20) 

P 

The Wilson loop (^|) is, naturally, preserved by this gauge transformation: it remains (— 1) 2J 
at p — > oo (for the representation labelled by spin J). 
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Having oriented the 'Higgs expectation value' Aq(oo) in one colour direction it is now 
possible to add up many vortices, each of them necessarily carrying a singular Dirac surface. 
Apart from the factor exp(— SVortex) the statistical weight of a vortex is determined by the 
fluctuation determinant, in particular by the zero modes of the solution. We expect three 
zero modes here, two of which are associated with the spatial position of the vortex center 
Zi, and one pure gauge cyclic mode associated with shifts in the 'time' direction. (The latter 
can be revealed if one uses a gauge with A explicitly dependent on time.) 

Denoting the (uncalculated) prefactor arising from the fluctuation determinant and from 
zero modes by K,(T,g 3 ) we can write the vortex partition function as that of a gas, 



^vortex = (/^ K eX P(~ ^vortex)) > ( 21 ) 



JV 

implying that at large T vortices are dilute and thus neglecting their interactions. It follows 



immediatelly from eq. (El]) that the spatial density of vortices (i.e. number per unit area) is 



n = Kexp(-S , vor t ex ). (22) 

At large T the numbers of noninteracting vortices inside any large area are Poisson-distributed. 
Therefore, the average of large Wilson loops can be calculated as 



lw\ X^\( U*A N " Are ^ N ^-n-Area J e - 2 "- Area at J = half-integer, 
<W0)-2_4(-l) j /yi 6 " 1 at J = integer, (23) 



N\ I 1 at J = integer, 

so that the string tension for half-integer representations is twice the vortex density [l[] , and 
zero for integer represenations. 

It should be recalled, however, that at T — > oo the 2 + 1 dimensional system reduces to 
two dimensions with a coupling constant g\ = g$T, and in two dimensions there is a trivial 
perturbative confinement with a nonzero string tension in any representation, 

a 2 = g 2 2 J(J+l). (24) 

The above vortex-induced string tension should be thus regarded as an exponentially small 
addition to eq. (p4|) : it leads to a deviation from the 'Casimir scaling' of eq. (pjj). It would be 
instructive to study the high-temperature 2+1 dimensional YM theory in lattice simulations, 
if only because one can learn to identify physical vortices by comparing objects found from 
maximal center gauge fixing with the continuum profiles presented here. Such experience 
might be useful in higher dimensions. 

As one lowers T the vortex action decreases; at T « g% the action becomes of the order 
of unity, therefore, the vortices are not suppressed anymore. Unfortunately, the quantitative 
theory fails at this point: first, because the effective action (H) is not accurate anymore, 
second, because the vortices become long and start to bend, third, because one can hardly 
neglect interactions between dense vortices. Much more serious efforts are needed to describe 
vortices in this case (if they exist). Nevertheless, to our mind it is useful to know that at 
least in the limit of high temperatures vortices do exist in a pure Yang-Mills theory, and 
their basic properties are established. 
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